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ABSTRACT 

The intervening large-scale structure distorts cosmic microwave background (CMB) anisotropics via 
gravitational lensing. The same large-scale structure, traced by dusty star-forming galaxies, also induces 
anisotropics in the far-infrared background (FIRB). We investigate the resulting inter-dependence of 
the FIRB and CMB with a halo model for the FIRB. In particular, we calculate the cross-correlation 
between the lensing potential and the FIRB. The lensing potential can be quadratically estimated from 
CMB temperature and/or polarization maps. We show that the cross-correlation can be measured with 
high signal-to-noise with data from the Planck Surveyor. We discuss how such a measurement can be 
used to understand the nature of FIRB sources and their relation to the distribution of dark matter. 
Subject headings: cosmology: theory - cosmology: observation - diffuse radiation - infrared:cosmology: 
weak lensing - cosmology: formation - galaxies: evolution 
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1. INTRODUCTION 

Dusty star-forming galaxies give rise to a far-infrared 
background (FIRB) (Puget et al. 1996; Fixsen et al. 1998; 
Dwek et al. 1998; Schlegel Finkbeiner & Davis 1998; La- 
gache et al. 1999; Guiderdoni et al. 1998; Blain 1999). 
Correlations in the large-scale structure traced by these 
contributing sources lead to correlated fluctuations in the 
FIRB (Bond et al. 1986; Scott & White 1999; Haiman & 
Knox 2000; hereafter HKOO; Knox et al. 2001; Maglioc- 
chetti 2001). At arcminute scales and more, fluctuation 
power associated with the source distribution can poten- 
tially be detected with Planck and other planned CMB ex- 
periments with channels at frequencies around and above 
300 GHz (Knox et al. 2001). 

The same large-scale structure that generates FIRB an- 
isotropy also generates anisotropy in the CMB in sev- 
eral ways. These include modifications due to scattering 
via free electrons in galaxy clusters, such as the thermal 
Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 1980), 
and modifications imposed by the time evolving gravi- 
tational field, such as the integrated Sachs- Wolfe effect 
(Sachs & Wolfe 1967). The large-scale structure mass field 
also deflects CMB photons propagating to us from the last 
scattering surface via the gravitational lensing effect (Sel- 
jak 1996). Since the lensing effect on CMB anisotropics is 
second order in temperature fluctuations, it induces non- 
Gaussian signatures in the temperature data; the cross- 
correlation between the lensing effect and other secondary 
anisotropics, such as SZ or ISW, contributes to the tem- 
perature bispectrum (Goldberg & Spergel 1998; Cooray & 
Hu 2000; Seljak & Zaldarriaga 1999). 

We extend previous discussions on correlations between 
the CMB and large-scale structure and consider the cross- 
correlation of CMB anisotropics and FIRB fluctuations. 
The FIRB contributes significantly at the high frequency 
end of certain CMB experiments, such as the 350 GHz, 
545 GHz and 850 GHz channels of the High Frequency In- 



strument (HFI) of the Planck Surveyor. Over this range 
of frequencies, and in regions of the sky with low galac- 
tic dust emission, the FIRB stands out as the dominant 
source of fluctuation power over a wide range of angu- 
lar scales (Knox et al. 2001). We model the expected 
cross-correlation between lensing potentials and the FIRB 
though a physically motivated semi-analytical approach 
involving the halo model (Scherrer & Bertschinger 1991; 
Seljak 2000; Scoccimarro et al. 2000; Cooray & Sheth 
2002). 

Though the lensing potential-FIRB correlation leads 
to a bispectrum or a three-point correlation function in 
CMB temperature and FIRB, we suggest a direct cross- 
correlation between lensing potentials and FIRB. Our sug- 
gestion follows from recent discussions in the literature on 
how to reconstruct lensing potentials associated with CMB 
lensing, especially through higher order statistics such as 
quadratic estimates optimized for the lensing extraction in 
temperature and polarization data (Zaldarriaga & Seljak 
1999; Hu 2001b; Hu & Okamoto 2001; Cooray & Kesden 
2002). We consider this possibility for planned missions 
such as Planck. 

Using reasonable assumptions, we show that Planck has 
sufficient sensitivity for a detection of the lensing-FIRB 
correlation. A detection of such a correlation would allow 
us to test how well the sources of the FIRB trace the large- 
scale mass distribution. While semi-analytic approaches 
such as the halo model suggest a high correlation, the ex- 
act correlation as measured will allow us to further refine 
details of these models and to understand certain physical 
properties of contributing FIRB sources. 

The layout of the paper is as follows. In Section 2 we 
present the halo model for FIRB fluctuations. In Sec- 
tion 3 we revisit the angular power spectrum of the FIRB 
and describe the cross correlation between the FIRB and 
a quadratic function of the CMB temperature map. The 
quadratic CMB statistic is an estimator for the lensing 
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potential. In section 4 we discuss associated errors and 
the extent to which the FIRB-lensing correlation can be 
detected. We conclude with a discussion of our results in 
Section 5. We refer the reader to Haiman & Knox (2000; 
hereafter HKOO) and Knox et al. (2001) for initial de- 
tails on our calculation of correlations in the FIRB. More 
details of the halo approach to large-scale structure are 
available in the recent review by Cooray & Sheth (2002). 
For a review of our observational knowledge of the FIRB 
see Hauser & Dwek (2001). While we provide a general 
derivation of the FIRB-lensing correlation, when illustrat- 
ing our results we assume a ACDM cosmological model 
with rirn = 0.35, flA = 0.65, = 0.05, h = 0.65. To de- 
scribe linear clustering, we use the transfer function given 
by Eisenstein & Hu (1999) and normalize fluctuations to 
COBE (Bunn & White 1997) such that erg = 0.86. 

2. HALO APPROACH TO THE FIRB 

As in HKOO, we write the FIRB Rayleigh-Jeans tem- 
perature (Trj = Ii,/{2kB){c/y)^) at frequency ;/ and in 
direction n as a sum of the mean and the fluctuation about 
the mean 

7Rj(n, v) = fnj{i^) + 5TRj(n, u) 

6j{rn, V, r) 



2A;bi/2 



dra{r)j{u, r) 



1 + 



(1) 

Here, r is the radial distance from us to a redshift of z and 
j{i',z) is the mean emissivity of the dust. Note that we 
have related the excess fluctuation in FIRB temperature 
to a fluctuation in the emissivity, 6j{rn, v, z). 

We adopt the HKOO model for the mean emissivity as 
a function of redshift and frequency. The key ingredients 
of this model are the history of ultra-violet radiation and 
dust production. Both of these are assumed to be pro- 
portional to a star-formation rate, which here and in the 
HKOO fiducial model we take to be that of Madau (1997). 
By further assuming the optical properties of the dust 
(Draine & Lee 1984) we then derive the comoving mean 
emissivity, jv{z). The dust and UV production propor- 
tionality constants are chosen so that the resulting FIRB 
agrees with inferences from COSiS/FIRAS data (Fixsen et 
al. 1998). 

To model fluctuations of the FIRB, we assume that 
the sources_ of this emissivity are galaxies so that 
5i{rn, V, z)/jiv, z) = (5gai where ngai = ngai(l + (5gai). Our 
modeling differs from that in HKOO in that we calculate 
the statistical relation of the galaxy number density field to 
the dark matter by making assumptions about how these 
galaxies populate dark matter halos; i.e., we use the 'halo 
approach' to large scale structure. 

The halo approach to large-scale structure involves a 
simplified description of the complex distribution of dark 
matter as consisting of a population of dark matter halos. 
The spatial distribution of any tracer field of the large- 
scale structure, such as galaxies, is described through 
the relation of the tracer field to that of the dark mat- 
ter halo distribution. The halo approach has now been 
widely discussed in the literature to understand the clus- 
tering of dark matter, galaxies and baryons (Scherrer & 
Bertschinger 1991; Scoccimarro et al. 2000; Cooray & 
Sheth 2002). The necessary inputs for the halo-based cal- 
culations come from either analytical tools, such as the 



dark matter halo mass function, or from numerical sim- 
ulations, such as the distribution of dark matter within 
each halo. 

The two-point correlation function of the object of inter- 
est in a halo model has two terms; the first is a two-point 
correlation function within each dark matter halo and the 
second is the two-point correlation function between points 
which are in different halos. The halo distribution is taken 
to trace the linear fiuctuations with biases which are pre- 
dictable from analytical arguments. The time evolution 
of the correlation function is naturally determined by this 
semi-analytic model since clustering evolution can be re- 
lated to the evolution of the halo mass function and any 
evolution of physical properties related to halos. As we 
find below, the halo model has the advantage that one 
does not need to specify a priori quantities such as the 
bias and evolution of fields that trace the mass distribu- 
tion of the large scale structure. 

Following Scherrer & Bertschinger (1991) and Scocci- 
marro et al. (2000) we can write the two-point function of 
the mass distribution as 

/0^(x - x', z) = 

j n{m, z)m^dm j (fyum{y, z)um{y -h x - x', z) 



X — Xi 



{mi,z)midmi J n{m2, z)m2dm2 J d^xiUmi C 

X I d.^X2Um^{x.' -X2,z)^{xi-X2-,mi.m2,z), 

where we have written out explicitly the redshift- 
dependence of the correlation function; hereafter, for sim- 
plicity, we usually drop this rcdshift-dcpendcncc when 
writing our expressions. In above, n{m) is the mass func- 
tion describing the number density of collapsed objects 
with given mass m: 



m?n{m) dm . , .dv 

= = i^fM — 

p m V 



(3) 



where p is the background density and j{v) is defined, 
following arguments given by Press & Schechter (PS; 1974) 
as 



Here, v is given by 



a'^(m, z) 

where S'^^{z) is the critical density for spherical cluster- 
ing and a-^{m, z) is the density fluctuation smoothed with 
a tophat fllter at a redshift of z corresponding to mass 
scale of TO. Useful fitting formula for 5sc{z) and related 
variables can be found in Henry (2000). The first term of 
equation (1) is the two-point correlation function within 
one halo, the second term is the two-point correlation func- 
tion between points being in two different halos and Uto(x) 
represents mass distribution in the collapsed region. 

We apply the NEW (Navarro Frcnk & White 1997) pro- 
file for the mass distribution in tlw^ collapsed region with 



Um{r) 



(6) 



(r/r«)(l + r/r,)2' 

where rg is the core radius and pg is the density at r^. 
We can also get from the known mean concentration 
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relation, c = r^ir/rs- We calculate the virial radius of each 
halo following spherical collapse arguments such that M = 
47rr^A(z)/9f,/3, where A(z) is the overdensity of collapse 
and ph is the background matter density today. We use 
comoving coordinates throughout. 

The mean concentration parameter c is given by (Bul- 
lock et al. 2001; Jing 2000) 

-I -0.13 

(7) 
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l + z 



m, (z) 



where m*(z) is the critical mass at given z. The density 
profile in Fourier space is 

2 sin kr p{r\m) 



u{k\m) 



drAnr 



(8) 



/g kr m ^ 

where r^ir is the virialized radius of the collapsed object. 

The Fourier transformation of equation (1) gives the 
power spectrum 

P^^{k) = J dmn{m) u{k\fn)'^ 

Pdmik) = / dmin(mi)^-u(k\mi) 
J P 

dm2n{m2)—r-u{k\m2) 
P 

X Phh(fc|"ii,m2). (9) 

Here, P/ift(fc|mi, TO2) is the halo-halo power spectrum 
which we take to be related to the power spectrum of the 
linear density field, P^(fc) by 

Phh(fc|mi,m2) = bi{mi)h^{m2)P^^(,k). (10) 
The halo biasing hi is (Mo et al. 1997) 

v-l 



bi{m) = 1 + 



(11) 



6sc{zi) 

We assume that the dust giving rise to the FIRB is all in 
galaxies. To describe clustering properties of these dusty 
galaxies, we modify the halo model for dark matter to 
describe point sources following previous discussions in the 
literature (Seljak 2000; Jing et al. 1998). 

An important ingredient is how these galaxies occupy 
dark matter halos. Following Scoccimarro et al. (2000); 
Cooray & Sheth (2002), we model this halo occupation 
number, the average number of galaxies in each dark mat- 
ter halo of mass m, through a simple relation of the form 

(A^gai(m)) = Agai {m/mof-^' . (12) 

When illustrating our results, we will take Agai ~ 
0.7, mo — 4.2 X 10^^ and Pgai = 0.8 consistent with what is 
suggested in the literature for the halo occupation number 
of blue, star-forming, galaxies (Scoccimarro et al. 2000; 
Cooray & Sheth 2002). Besides these parameters, we also 
set a lower bound on the halo mass, at rric = 10^^ M^/h, 
in which dusty galaxies may exist. This accounts for the 
fact that baryons may not collapse and cool to form stars 
in low mass halos as well as to account for the fact that 
star-formation may not be efficient due to processes such 
as supernova winds which may expel any remaining gas in 
such small mass halos (Benson et al. 2001; Kauffmann et 
al. 1993). 




Fig. 1. — Galaxy bias as a function of k for our fiducial halo 
model at 2 = 0,1,2 (upper panel) and galaxy-dark matter density 
contrast correlation coefficient (lower panel). 



Note that we have described the halo occupation number 
of sources that contribute to FIRB with a description that 
is consistent for blue galaxies at low redshifts. While these 
late-type galaxies are expected to trace dusty-starbursts 
that form the FIRB, our choice of parameters is also con- 
sistent with another consideration. The non-linear cluster- 
ing of sources in the IRAS Point Source Redshift Survey 
Catalog (PSCz), based on 60 /zm fluxes, is well produced 
by a halo occupation number consistent with the above set 
of parameters for p, Agai, and ttiq (Scoccimarro et al. 2000; 
Cooray 2002). These sources are also likely to trace the 
dusty-starbursts that contribute to FIRB. An additional 
complication, which we return to later, is if FIRB sources 
are associated with either mergers or active galactic nuclei. 

To describe galaxy clustering at the two-point level we 
also need the second moment, (-/Vgai(A'gai— 1)) of the galaxy 
halo occupation distribution. Given our limited knowl- 
edge, we take the simplest approach involving a Poisson 
distribution such that (iVgai(A^gai - 1)) = {Ngai{m))2. This 
assumption is consistent with numerical simulations espe- 
cially at the high mass end while at the low mass end of 
the mass function, when {Ngs,i{m)) < 1, evidence exists for 
departures from a Poisson distribution (Scoccimarro et al. 
2000; Cooray & Sheth 2002). 

With information related to the halo occupation num- 
ber, we can write the power spectrum of dusty galaxies, 
again with two terms: 



P^tiik) 



P^hik) 



dmn{m 



/(iVgalM^ 
V %al(™) 



u{k\mY 



dmin(mi)^ — —u(k\mi) 

"gal (mi) 
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J dm2n{m2) 



(iVgal(m2)) 



u{k\m2) 



"gal(m2) 

X Phh{k\mi,m2). (13) 
Here, p accounts for the possibility that at least one galaxy 
may reside in the center of each halo. Thus, we set p = 1 

when (iVgai) < 1 and otherwise p = 2. Here, figai is the 
mean density of FIRB sources given by 



"gal = / dmn{m){Ngi,i{m)) 



(14) 



The angular power spectrum of FIRB fluctuations can now 
be obtained by projecting the three-dimensional dusty 
galaxy power spectrum along the line of sight following 
standard techniques such as the Limber approximation. 
The cross correlation between lensing and the FIRB, how- 
ever, depends on the correlation of dusty galaxy and dark 
matter distributions. The halo model provides a simple 
calculation of this cross-correlation and we write 

Pdm-gal{k) = Pdm-gali^) ' ^ dm-ga 

,m(A^gai(m)) 

-gal V 



plh 



dm— gal 



J P ngs,i[m) 

(k) — / dmin{mi)^^u(k\mi) 
J P 

X / dm2n{m2)- 



\ni2) 



(15) 



"gal(m2) 

X Phh{k\mi,m2). 
The halo approach can be used to estimate the dusty 
galaxy bias, &gai, and correlation, rgai, both as a function 
of scale, k, and redshift z. We define these two quantities 
as 



bga.\{k,z) = 



Pdmik, Z) 



Pd: 



Im—gal 



(16) 

V-Pgal-Pgal 

and plot them in Fig.l as a function of k for three redshifts 
(2; = 0, 1 and 2). Note that &gai(^) is constant at the large 

spatial scales corresponding to the linear regime. This 
large-scale bias, under the halo model, is simply 

(A^gai(m)) 



bga.\{k) — / dmn{m)bi{m) 



(17) 



At non-linear scales bias increases due to the choice of 
p = 1 when (iVgai) < 1. The behavior of rgai(fc) can also be 
understood from the halo model. At linear scales, galaxies 
trace the dark matter such that they perfectly correlate 
with each other. At non-linear scales, rgai(fc) is greater 
than unity since the galaxy power spectrum is defined such 
that, following the standard definitions, it does not include 
the shot noise part associated with the finite number of 
galaxies (Seljak 2000). 

3. ANGULAR POWER SPECTRA 

In this paper we arc interested in measuring the cross- 
correlation between FIRB fluctuations and the lensing ef- 
fect on CMB anisotropics. For the purpose of this dis- 
cussion, we will describe lensing through the associated 
projected potential, (j), along the line of sight 



(f){n) = -2 dr 



dA {rp - r) 
dA{r)dA{ro) 



$(r, mr) , 



(18) 




2000 
r (h 



4000 
•Mpc) 



6000 



Fig. 2. — Integrand of equation (27). Here wc 
plot the integrand of equation (27) broken up into 
two factors, W''"{k,u, z)/r,y P^^_gai{l/r) (light) and 

(k, u, z)/ry/ P^^_galil/^) (hsa'Vy). Top panel shows I = 100 
(dashed), I = 400 (solid) and I = 1000 (dotted) at u = 545 GHz 
and bottom panel shows u = 353 GHz (dashed), fc" = 545 GHz 
(solid) and u = 852 GHz (dotted) at I = 400. 



where ro = r{z = 1100) corresponds to the radial distance 
to background CMB fluctuations which act as a source for 
lensing and <I> is the gravitational potential field. Here, 
dA is the comoving angular diameter distance and, in a 
spatially-flat cosmology, dA ^ r. Note that the deflection 
angle of a CMB photon is given by the gradient of this 
projected potential, a{n) = W(f>(h). 

First, we will consider the three power spectra associ- 
ated with the projected lensing potential, FIRB fluctua- 
tions at frequency v, F^, and the cross-correlation between 
the two. We use flat-sky coordinates throughout the paper 
and write the Fourier transform of quantities as 

X{\) = J d^eX{e)e-'^-^ , (19) 

where X = F,(j) for FIRB and lensing potentials respec- 
tively. We define power spectra of these quantities and the 
cross-power spectrum such that 

{FAI)F,{1')) ^ {2nrSn{\ + l')Cf^{i^), 

mW)) = {2nfSn{l + l')Cf1', 
{F^mr)) = {27rf6n{l + l')Cf^{>^). 



(20) 
(21) 

(22) 
(23) 

Wc can write Fourier coefficients for fluctuations in either 
the FIRB or lensing potential as 

X(l)= J^Jd^eJ drW^^(A,r)<5x(fc,r)e-'('-'''-)■^ 

(24) 

Following Knox et al. (2001), we write the FIRB weighting 
function at a given frequency ly as 

W^{z; v) = a{z)j{iy, z)/{2kB){c/yf and 

5f = 5gal. (25) 
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Similarly, for lensing potentials. 



100 



Ho 

k 



dAiro - r) 



and 



a{r) \ k J dA(rjdA{ro) 

^0 = (^dm- (26) 

We have used Poisson's equation to replace the potential 
field with the dark matter density field, resulting in the 
fc-dependence of the window function. 

As written in equations (25) and (26), and discussed ear- 
lier, FIRB fluctuations trace that of contributing galaxies, 
while lensing traces fluctuations in the dark matter field. 
For illustration purposes, we plot both FIRB and lensing 
weight functions in Fig. 2 as a function of multipole, or 
angular scale, and frequency. Note that the lensing weight 
function, while peaking at nearly half the angular diameter 
distance to background source, is a broad function in red- 
shift space. This allows lensing potentials to be correlated 
with a wide variety of large-scale structure tracer fields. 
The main result from this comparison is that the lensing 
and FIRB weight functions overlap with each other over 
a wide range in radial distance space. As we find later, 
this substantial overlap results in a significant correlation 
between lensing potentials and FIRB fluctuations. 

Using the weighting functions in radial space, we can 
now calculate the angular power spectra of FIRB and 
lensing potential fluctuations as well as the cross angular 
power spectrum between the FIRB and lensing potentials. 
We make use of the the Limber approximation (Limber 
1954; Kaiser 1992) and write 



Cf 



dr 

"-A 
dr 

< 



[Ty'^"(fc,z)]^Pd. 



I 

dA 



dA ' 



m— gal 



dA ' 



(27) 

for lensing-lensing, FIRB-FIRB and FIRB-lensing power 
spectra. 

In Fig. 3 we show the angular power spectrum of the 
FIRB for our fiducial model at 545 GHz. Here, we have 
calculated five predictions for the FIRB signal including 
two constant bias models that trace the non-linear dark 
matter power spectrum under the formulation given by 
Peacock & Dodds (1996) with bias equal to 1 (lower curve) 
and 3 (upper curve). The solid curve in the middle is the 
prediction for the FIRB with our halo model. Because 
u = 545 GHz is much greater than the frequency of peak 
intensity for a 2.73 K black body, the CMB signal is sig- 
nificantly smaller than the FIRB signal. 

For comparison, in Fig. 3, we show the two contributions 
to the FIRB power spectrum under the halo model involv- 
ing 1- and 2-halo terms. As shown, the FIRB clustering at 
tens of arcminute scales and above is described by correla- 
tions of FIRB sources in different halos while at arcminute 
scales, the 1-halo term due to correlations of sources within 
the same halo dominates. At these small scales, however, 
instrumental noise starts to contribute rapidly contami- 
nating the FIRB clustering at smallest scales. 

In Fig. 4, we show the cross-correlation between FIRB 
fluctuations at 545 GHz and the lensing potential. The 
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Fig. 3. — Angular power spectra at 545 GHz of our fiducial FIRB 
model (solid with error boxes), CMB (dot-dashed), galactic dust 
averaged over cleanest 10% of the sky (short/long dash) and Planck 
noise (dashed). The curves labeled 1-h and 2-h show the 1- and 
2-halo contributions to the FIRB power spectrum. Error boxes are 
for measurement of FIRB power spectrum by Planck using only 
the cleanest 10% of the sky and including the contributions from 
unsubtracted galactic dust. The dashed lines below the fiducial 
FIRB model are for angular power spectra for FIRB models with 
rrxc — > IO^^Mq/Zi (lower) and Pg^i 0.9 (upper). Dotted lines are 
constant-bias FIRB models with 6 = 1 (lower) and b = 3 (upper). 



curves shown here follow those of Fig. 3. As shown, cross- 
correlation in the multipole range of interest is dominated 
by the 2-halo term. This captures the large scale linear 
correlations between FIRB sources and dark matter in dif- 
ferent halos instead of the clustering of FIRB sources and 
dark matter in the same halo. This dependence on the 
two-halo term is understandable from the weight functions 
plotted in Fig. 1. As shown there, the cross-correlation 
peaks where the weight functions overlap the most and 
this is at radial distance corresponding to redshifts of 1 
and higher. At such distances, the non-linearities, cap- 
tured by the 1-halo term, do not dominate. 

The strength of the correlation between and F is quan- 
tified by the cross-correlation coefficient 

iF0 



.(0 



(28) 



In Fig. 5 we show that the correlation is moderate, as 
expected from the rough similarity of the weighting func- 
tions. Note that if fraction 1 — e of the FIRB signal is 
completely dependent on the lensing signal (i.e., can be 
predicted from it) then we expect 1 — rcorr — 0.5e^. A 
significant uncorrelated component still results in a rcorr 
near unity; e.g., e = 0.4 — > rcorr = 0.92. 

In addition to cross-correlation, we can also define a bias 
factor in Fourier space relating FIRB fluctuations and dark 
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Fig. 4. — Angular power spectrum of FIRB-CMB lensing cross- 
correlation function at 545 GHz for same FIRB models as in Fig. 3. 
The curves labeled 1-h and 2-li show the 1- and 2-halo contributions 
to the cross power spectrum. Note that in the range of angular scales 
of interest to Planck, all contributions effectively come from the 2- 
halo term or the large scale correlations between FIRB sources and 
dark matter between different halos. Error boxes assume use of 10% 
of the sky, include the contributions from unsubtracted dust and do 
not include polarization information. Solid (dashed) error boxes are 
for Planck (no) noise. 
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Fig. 5. — Top: The correlation coefficient (see equation (28)) at 
545 GHz. (a) Error boxes assume the observation of the cleanest 
10% of the sky with Planck noise and include the contributions 
from unsubtracted dust, (c) For comparison in this panel we show 
the error boxes obtained when thereis no dust contribution to the 
FIRB map (solid line) and, in addition, no noise contribution to the 
lensing reconstruction such that Nf''' = (dashed line). Bottom: 
The bias coefficient for FIRB fluctuations relative to that of the 
lensing potential field. The error boxes follow those in the top panel. 



matter traced by lensing potentials, such that 
b{l)-- 




(29) 



We also plot this bias factor and its associated errors on 
the bottom panel of Fig. 5. 

4. ERROR FORECASTS 

Here we consider expected measurement errors on the 
various signals predicted in the previous sections and 
shown in Figures 3, 4, and 5. First, we consider the ex- 
pected errors on the detection of lensing potential power 
spectrum and the FIRB power spectrum. Then we will 
discuss the associated errors on the determination of the 
cross-correlation between FIRB and lensing. 

4.1. Lensing 

In the case of lensing power spectrum, we assume that 
the deflection angle associated with the CMB lensing can 
be constructed from quadratic statistics involving temper- 
ature fluctuations (Hu & Okamoto 2001; Cooray & Kesden 
2002). 

To understand the mechanism involved, first notice 
that lensing involves a remapping of the CMB tempera- 
ture (Seljak 1996; Zaldarriaga 2000; Hu & Okamoto 2001; 
Cooray & Kesden 2002) such that 
e(n) = e[h + \/(j){h)] 

= e(n) + V,0(n)V*e(n) 

+ iv,0(n)Vj</)(n)V*V^e(n) + . . . (30) 



where 0{n) is the unlensed primary component of CMB, 
the contribution at the last scattering surface, and 6(n) is 
the total contribution observed today after modifications 
due to gravitational lensing deflections during the transit. 
To extract information related to lensing potentials, now 
consider a simple quadratic statistic involving the squared 
temperature map, which in Fourier space can be written 
as 

e'(i) = / - ii) , (31) 

where ly is a filter that acts on the CMB temperature and 
can be derived by the requirment that an ensemble aver- 
age of this quadratic statistic returns an estimator for (f), 
(0^(1)) = (f>{l), with a minimum noise contribution. Tak- 
ing the Fourier transform of equation (30), and substitut- 
ing for 0(1) in equation (31), one can obtain the required 
filter as 

^^[MiCf -1-1- (l-l^Ci? , ,] 

^•--ii '-'li-iil 
with a normalization given by 

dH, [MiC« + l-(l-li)Ce_,^|]^ 



(2^)2 



2^tot ^tot 
h |1— li 



(33) 



Here, Cf is the power spectrum of unlensed primordial 
temperature fluctuations, as measured at the last scat- 
tering surface, while C*°* acounts for all contributions 
to the temperature power spectrum such that C*°* = 



lonscd-CMB 



Here, 



lensed-CMB 



is the 
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CMB power spectrum after gravitational lensing, C^^o'se jg 
detector and instrumental noise contributions, and Cf"'^'^ 
describes the contribution due to any foregrounds and sec- 
ondary effects, such as the SZ effect. 

Note that the variance of the filtered quadratic temper- 
ature map. 



cf 



■Nf 



(34) 



returns the lensing potential power spectrum plus an as- 
sociated Gaussian noise bias given by Nf'^. We note that 
even in the case of a CMB experiment with no detec- 
tor noise, the reconstructed potential, under the quadratic 
statistic techhnique described above, is biased as Nf^ has 
contributions from primordial CMB fluctuations. This im- 
plies that the correlation coefficient defined in equation 
(28) is not the correlation coefficient one measures between 
the FIRB map and the reconstructed map, as described 
above, even for an ideal experiment. This noise bias, how- 
ever, can be eliminated possibly through a different re- 
construction technique involving an unbiased estimator of 
the lensing deflection potential. In order to generalize our 
discussion, we ignore this noise-bias and treat Nf"^ as the 

error contribution associated with the reconstructed Cf'^. 
Thus, correlation coefficients we plot in Fig. 5 are those 
that one would measured by correlating the FIRB map 
with an unbiased estimator of the deflection angle. 

In addition to the case with temperature we have just 
discussed, the polarization field can also be used to recon- 
struct (j) (Guzik, Seljak, & Zaldarriaga 2000; Benabed et 
al. 2000; Hu & Okamoto 200f ; Cooray & Kesden 2002), 
though we perform no error forecasts for polarization data 
here. 

When estimating iV^'^'^, we consider two experiments cor- 
responding to Planck and a no-detector noise, cosmic vari- 
ance limited, experiment. In the case of Planck, wc use 
tabulated values for its detector sensitivities and resolu- 
tions, or beam sizes, following Cooray et al. (2000) to cal- 



culate C," 



In the case of the cosmic variance limited 



experiment, we set C" 



and take a maximum I of 



3000. to which we integrate equation (33). In both these 
experimental possibilities, we set Cf"^^ = 0. 

4.2. FIRB 

On the FIRB side, we expect a significant contribution 
to the angular power spectrum at 545 GHz from galactic 
dust. The level of the dust power spectrum shown in Fig. 3 
is for that in the cleanest 10% of 5° x 5° patches (Knox et 
al. 2001). The average level in the cleanest 50% is about 
10 times higher. 

We approximate the effect of dust contamination by in- 
cluding it as an iinsubtracted noise source with a power 
spectrum given by Cf'. The error on the FIRB-FIRB 
power spectrum is 

1 



AC, 



FF 



^{21 + l)/,ky 



FF.tot 



where following Knox et al. (2001) 



cj^''°'={cr +ci 



■N, 



FF\ 



and 



I(l+l)a2/81n2 



(35) 



(36) 



Here, Nf^ is due to instrumental noise in the beam- 
deconvolved maps. For example, the Planck 545 GHz 
channel has a = 5.0' and (on average) Atrj = 
13/iK — arcmin. Throughout we use the u = 545 GHz 
channel as our probe of the FIRB though valuable, and 
complementary, information does lie in other channels at 
both higher and lower frequencies. 

Equation (35) is formally correct for a single-frequency 
measurement if the power spectrum of the dust is perfectly 
known and the dust is a statistically isotropic Gaussian 
random field. These conditions do not hold; this equa- 
tion only serves as a rough guide, with dust increasing 
the errors when Cf > Cf^ and having little effect when 
cf « Cf^. Throughout we restrict our forecasting to 
results from the cleanest 10% of the sky, using the appro- 
priate dust level. 

4.3. Lensing-FIRB 

Following our discussion on lensing reconstruction, one 
can use the estimator for the deflection angle, or the 

projected potential, to measure the cross-correlation with 
FIRB fluctuations. Thus, similar to equation (34), we can 

write (e2(l)F(l')) = (27r)2^D(l + l')Cf as an estimate of 

the cross-power spectrum between lensing potentials and 

FIRB fluctuations. In Fig. 4 we plot Cf**. Note that 

while Cf"^ is biased under quadratic temperature recon- 
struction techniques, the cross-correlation power spectrum 
between lensing potential and FIRB is unbiased; this is 
due to the fact that FIRB fluctuations are assumed to 
be uncorrelated with unlensed temperature anisotropies 

(e(i)ni'.)> = 0. 

Following our discussion of errors on lensing and FIRB 
fluctuation power spectra, we can write the expected error 
on Cf"^ as 

1/2 



1 



{21 + 



+ (Cf^'*°*)(C"''°') 



FF.tot-, 



where 



,tot 



(37) 

is the total contribution in 

the <f) power spectrum. 

The associated error on the correlation coefficient and 
the bias factor can be calculated by propagating errors 
related to Cf'^, Cf^, and, Cf'^. For reference, we repro- 
duce the flnal result for these errors here. In the case of 




while, in the case of 6(0 yJcf^/Cf'^^ : 



6^(0 



2(2/ + l)Uy [ cf 



A.tot 



c, 



FF.tot 



^FF 



(39) 
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5. DISCUSSION 

We now discuss the main results of our paper. We have 
introduced a scini-analytic approach to describe the FIRB 
fluctuations. The model involves a distribution of dark 
matter halos populated by dusty, starforming galaxies. 
The basic ingredients of the calculation include properties 
of this dark matter distribution, such as the mass func- 
tion and clustering properties of dark matter halos, and 
properties of the sources responsible for the FIRB. To de- 
scribe these sources we have introduced a relation, the halo 
occupation number, that attempts to capture how FIRB 
sources populate dark matter halos. These model can be 
used to calculate measurable properties of the sources such 
as their bias as a function of redshift and their spatial cor- 
relation with the underlying dark matter as a function of 
redshift. 

The approach presented here differs from previous at- 
tempts to model the fluctuations in the FIRB. HKOO 
used several biasing models, all of which were scale- 
independent. The one most similar to our calculation here 
assumed all the FIRB somces were in 10^^ Mq halos and 
used the biasing prescription of Mo & White (1996). Be- 
cause we include a range of halo masses extending below 
10"'^^ Mq and because we take the FIRB mean to be \/2 
lower than HKOO did our power spectra have smaller am- 
plitudes. The fluctuation power we predict is also smaller 
than that in Knox et al. (2001) who used a constant 6 = 3 
amplification of the non-linear power spectrum calculated 
using the fitting formula of Peacock & Dodds (1996). 

Even though we predict a smaller power spectrum than 
in the constant bias models we still expect FIRB fluctua- 
tions to be detectable in upcoming high frequency experi- 
ments such as Planck. Any multifrequency detection will 
allow the testing of our models and hopefully the extrac- 
tion of interesting properties of our universe such as the 
evolution of the star-formation rate. We do not discuss 
such possibilities in detail here as they have been already 
addressed in Knox et al. (2001) and references therein. 

In this paper we have gone beyond tlic^ FIRB fluc;tuation 
power and focused on the cross-correlation between FIRB 
sources and the dark matter as traced by the intervening 
lensing potential that deflected CMB photons propagating 
to us. Note that previous studies have already suggested 
that experiments such as Planck will detect the lensing po- 
tential power spectrum with considerable signal-to-noise 
ratios (Zaldarriaga & Seljak 1999; Zaldarriaga 2000; Hu & 
Okamoto 2001; Cooray & Kesden 2002). In addition to the 
separate detections of FIRB and projected lensing poten- 
tial power we have suggested that one can also perform a 
combined study to measure the cross-correlation between 
FIRB fluctuations and the lensing potential. 

We have discussed an extension of the estimators sug- 
gested for lensing reconstruction in CMB data. Using 
quadratic statistics one can deflne an estimator of the de- 
flection angle, which can then be directly correlated with a 
high frequency map where FIRB fluctuations are expected 
to dominate. Such a direct approach avoids complications 
associated with other methods for extracting the lensing- 
FIRB cross power spectrum such as the direct measure 
of the three point function. As discussed in the literature, 
the FIRB-lensing correlation leads to a three-point correla- 
tion function, or a bispectrum in Fourier space (Goldberg 



& Spergel 1998; Cooray & Hu 2000). Note that measur- 
ing the full configuration dependence of this bispectrum 
is difficult and currently limited by computational meth- 
ods and measurement techniques. While improvements 
are expected, our suggested technique avoids these issues 
by combing the different three point function configura- 
tions in a particular way. 

As shown in Fig. 4, we find a considerable correlation 
between FIRB fluctuations and the lensing potential. As 
illustrated in Fig. 2, this is due to the broad overlap be- 
tween the respective radial weighting functions. Such a 
high correlation coefficient also leads to the conclusion 
that upcoming CMB experiments will be able to detect 
the cross power spectrum between FIRB and lensing with 
high signal-to-noise. For Planck, we expect the cumula- 
tive signal-to-noise to be of order 40 in the p = 545 GHz 
channel. 

The high signal to noise expected for the detection of 
the cross correlation and of the FIRB and lensing power 
spectra suggests that we will be able to test our under- 
lying model for the clustering of FIRB sources. As we 
discussed our approach needed two main ingredients, the 
bias of FIRB contributing sources which was calculated 
in the halo model, and the star formation rate which is 
needed to obtain FIRB weighting function. In turn, the 
halo approach depends on how starforming dusty galaxies 
populate dark matter halos. As shown in Fig. 3 and Fig. 
4 we expect that diff'erent values of the model parameters 
can be distinguished with the level of noise expected for 
Planck. 

In Fig. 5, we summarize our results with respect to how 
well the projected correlation coefficient, rcorr(Oi ^^^^ 
bias factor, b{l), can be measured with upcoming Planck 
data. While Planck allows a reliable detection of the cor- 
relation coefficient as well as bias out to a multipole of 
^ 1000, there are further improvements one can hope 
to achieve in the post-Planck era. As summarized in 
Fig. 5(c), in the case of the correlation coefficient, the error 
is dominated by the uncertainty associated with the lens- 
ing reconstruction. With improved data involving better 
angular resolution and noise one can hope to reach the lim- 
iting case of a perfect experiment which we have demon- 
strated with a dashed line Fig. 5(c). 

A few words of caution are required at this point. Note 
that we have described the FIRB and its fluctuations as 
coming only from dusty starburst galaxies. In addition 
to the stellar contribution, there is also some from dust 
associated with torii surrounding active galactic nuclei. 
While the fraction of AGN contribution at far-infrared 
wavelengths could be as high as 30% to 40% (Almaini 
et al. 2001; Risaliti et al. 2002), we have limited informa- 
tion on their redshift evolution so we have not included 
them in our calculation. The inclusion of AGNs within 
the halo model is trivial, the problems arise when trying 
to calculate the weighting function. 

Another caveat comes from our simple description of 
how submm sources populate halos. While we have de- 
scribed them through a halo occupation number, the 
physics is likely to be more complicated especially if dusty 
star formation is associated with mergers. With a merger 
rate proportional to N^i^i{m), the effective occupation 
number will probably depend on a higher power of mass 
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than considered here. In such a scenario, we expect source 
bias to be larger with an increase in clustering power com- 
pared to the results presented here. Such an increase 
should be detectable and constrained using observational 
data. 

While these issues may complicate the direct interpreta- 
tion of the observations, a reliable detection of the FIRB- 
lensing correlation can allow one to introduce more so- 
phisticated analysis techniques to try to understand these 
complicating factors. For example, the two power spectra 
and the cross-power spectrum can be combined and writ- 
ten as = b2{l)Cf* and Cf^ = b{l)r{l)Cf^, which can 
be inverted with appropriate techniques to obtain r{k, z) 
and h{k, z). Such an inversion requires accurate informa- 
tion on the FIRB radial weight function which can be ob- 
tained observationally through the redshift distribution of 
sources that contribute to FIRB. While unresolved surveys 
such as Planck will not allow such studies, in the future, 
targeted studies of resolved FIRB sources, over small but 
representative patches of sky may provide the necessary 
information. Constraints on this emissivity as a function 
of redshift from current data are discussed in Gispert et al. 
(2000), Chary & Elbaz (2001) and Takeuchi et al. (2001). 
With an inversion of two-dimensional clustering to three- 
dimensions one can eventually constrain various aspects 
of the halo model such as the occupation number (Cooray 
2002). While previous studies have motivated the use of 
FIRB fluctuations alone to understand an important as- 
pect of the large scale structure and its evolution history, 
we also suggest that cross-clustering aspects, such as be- 
tween FIRB and dark matter as traced by gravitational 
lensing of CMB also plays an important role. 
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Grant NAG5-11098. MZ is supported by NSF grants AST 
0098606 and PHY 0116590, and by the David and Lucille 
Packard Foundation Fellowship for Science and Engineer- 
ing. This research was supported in part by the NSF under 
Grant No. PHY99-07949. 
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